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Neutral pion production in central 158-A GeV 208 Pb+ 208 Pb collisions has been studied in the 
WA98 experiment at the CERN SPS. The ty° transverse mass spectrum has been analyzed in terms of 
a thermal model with hydrodynamic expansion. The high accuracy and large kinematic coverage of 
the measurement allow to limit previously noted ambiguities in the extracted freeze-out parameters. 
The results are shown to be sensitive to the shape of the velocity distribution at freeze-out. 

25.75.Dw, 24. 10. Pa 



1 



Heavy ion reactions at sufficiently high energies pro- 
duce dense matter which may provide the necessary con- 
ditions for the transition from a hadronic state to a de- 
confined phase, the Quark-Gluon Plasma. Since a finite 
thcrmalizcd system without external containment pres- 
sure will necessarily expand, part of the thermal excita- 
tion energy will be converted into collective motion which 
will be reflected in the momentum spectra of the final 
hadrons. The dynamics of the expansion may depend on 
the presence or absence of a plasma phase. The strongly 
interacting hadrons are expected to decouple in the late 
stages of the collision. Their transverse momentum spec- 
tra should therefore provide information about the con- 
ditions of the system at freeze-out, in particular about 
the temperature and collective velocity of the system, if 
the thermal assumption is valid. 

The application of a thermal description is non-trivial. 
There is no reason to believe neither that chemical and 
kinetic freeze-out should be identical, nor that there 
should be unique thermal freeze-out temperatures for all 
hadrons, nor unique chemical freeze-out temperatures for 
all flavour changing reactions. It is likely that chemical 
equilibrium is not fully attained (see e.g. |Q), implying 
that chemical parameters will also influence momentum 
spectra through contributions from decays of heavier res- 
onances. Furthermore, it is not obvious that this problem 
should have a stationary solution since particle emission 
will occur throughout the full time evolution of the colli- 
sion and so, in principle, would require a full space-time 
integration with varying parameters. 

Most attempts to extract freeze-out parameters from 
experiment assume local thermal equilibrium and fit pa- 
rameterizations of hydrodynamical models to the exper- 
imental distributions Already the earliest analy- 
ses Q| noted ambiguities in fitting the hadron transverse 
mass spectra due to an anti-correlation between the fitted 
temperature, T, and transverse flow velocity, fix- 

Two-particle interferometric (HBT) measurements 
provide information on the spatial and temporal extent 
of the emission volume, but are also sensitive to the col- 
lective motion of the source (see e.g. Within a 
hydrodynamical parameterization of the source at freeze- 
out, the transverse two-particle correlations have been 
shown to be sensitive only to the ratio f3 T /T S. Hence 
HBT analyses have a /3t —T ambiguity which is roughly 
orthogonal to that resulting from fits to the single par- 
ticle spectra. This fact has recently been used by the 
NA49 collaboration to constrain the freeze-out param- 
eters to lie within the region (/3t) — 0.55 ± 0.12 and 
T = 120 ± 12 MeV for central Pb+Pb collisions §. Al- 
ternatively, a recent analysis of tt + ,K + , and K~ distri- 
butions and 7r + and 7r~ two-particle correlations mea- 
sured by the NA44 collaboration for central Pb+Pb col- 
lisions using a 9-parameter hydrodynamical model fit JTcj ] 
gave freeze-out parameters of (fir) — 0.443 ± 0.023 and 
T = 95.8 ± 3.5 MeV. These analyses suggest that a sin- 



gle set of freeze-out parameters can describe the hadron 
single particle distributions and two-particle correlations, 
with moderate temperature and large transverse flow ve- 
locity. 

On the other hand, various thermal model analyses 
of particle production ratios, especially strangeness pro- 
duction (see e.g. Ref. jllj for a recent summary), have 
indicated rather high chemical freeze-out temperatures. 
Use of integrated yields in these analyses allows to obtain 
conclusions on the temperature which are insensitive to 
the amount of transverse flow. In a recent analysis of re- 
sults at SPS energies, including Pb+Pb collisions, good 
agreement is obtained if partial strangeness saturation is 
assumed with a chemical freeze-out temperature of about 
180 MeV |J]. 

A successful thermal interpretation of relativistic 
heavy ion collisions must provide an accurate descrip- 
tion of the pion spectra since pions provide the "thermal 
bath" of the late stages the collision. In this letter we 
discuss the extraction of thermal freeze-out parameters 
from the neutral pion transverse mass distribution for 
central 158-A GeV 208 Pb+ 208 Pb collisions. These data 
provide important constraints due to their accuracy and 
coverage in transverse mass. The analysis of the ir° spec- 
trum, within a particular hydrodynamical model, reveals 
the importance of the shape of the velocity distribution 
at freeze-out. The default shape, derived from a Gaus- 
sian spatial distribution, favors a large thermal freeze- 
out temperature, similar to temperatures extracted for 
chemical freeze-out, but in contradiction to conclusions 
obtained based on analyses of limited coverage particle 
spectra and HBT results 

The CERN experiment WA98 consists of large 

acceptance photon and hadron spectrometers together 
with several other large acceptance devices which allow 
to measure various global variables on an event-by-event 
basis. The results presented here were obtained from 
an analysis of the data taken with Pb beams in 1995 
and 1996. The 10% most central reactions (cr centra ; ks 
630 mb) have been selected using the transverse energy 
Et measured in the MIRAC calorimeter. 

Neutral pions are reconstructed via their 77 de- 
cay branch using the WA98 lead-glass photon detector, 
LEDA, which consisted of 10,080 individual modules with 
photomultiplier readout. The detector was located at 
a distance of 21.5 m from the target and covered the 
pseudorapidity interval 2.35 < r\ < 2.95. The general 
analysis procedure, described in Jr|, is similar to that 
used in the WA80 experiment 0. The momentum dis- 
tributions are fully corrected for geometrical acceptance 
and reconstruction efficiency. The systematic error on 
the absolute yield is « 10% and increases sharply below 
Pt = 0.4GeV/c. An additional systematic error origi- 
nates from the uncertainty on the momentum scale of 1%. 
The influence of this rises slowly for large px and leads 
to an uncertainty on the yield of 15% at pr — 4 GeV/c. 
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FIG. 1. Transverse mass spectra of neutral pions in cen- 
tral collisions (10% of min.bias cross section) of 158 AGeV 
Pb+Pb. The invariant cross section of neutral pions is com- 
pared to a fit using a hydrodynamical model including 
transverse flow and resonance decays, with the direct produc- 
tion and the contribution of p decays and all other resonances 
shown separately. The ratio of the fit to the data is shown in 
the inset, mo is the 7r° mass. 

The measured neutral pion cross section from central 
Pb+Pb reactions as a function of tut — mo is shown in 
Fig. 0. Included is a fit with a hydrodynamical model 
H including transverse flow and resonance decays. This 
computer program calculates the direct production and 
the contributions from the most important resonances 
having two- or three-body decays including pions (p, Kg, 
K*, A, £ + A, rj, u>, r/). The code, originally intended 
for charged pions, has been adapted to predict neutral 
pion production. The model uses a gaussian transverse 
spatial density profile truncated at 4<r. The transverse 
flow rapidity is assumed to be a linear function of the 
radius. For all results presented here, a baryonic chemical 
potential of hb = 200 MeV has been used. The results 
are not very sensitive, however, to the choice of \ib for 
the tut — mo region considered here. 

This model provides an excellent description of the 
neutral pion spectra with a temperature T — 185 MeV 
and an average flow velocity of (/3t) = 0.213. These 
values are very similar to the parameters obtained with 
similar fits to neutral pion spectra in central reactions of 
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FIG. 2. The local inverse slope of the transverse mass 
spectrum of neutral pions in central collisions of 158 AGeV 
Pb+Pb. The measured results (solid points) are compared 
to the hydrodynamical model best fit result (solid line; 
T = 185 MeV and (fir) = 0.213, also shown in Fig. 1) and to 
the other results given in table H. 

32 S+Au 0. The 2cr lower limitQ on the temperature is 
rpiow _ yj-^ MeV and the corresponding upper limit on 
the flow velocity is {(3j? p ) = 0.253. 

The observed curvature at low bit is largely a result of 
resonance decay contributions. Performing a fit with only 
the direct contribution leads to T = 142 MeV and (/3t) = 
0.301, with corresponding 2er limits of T low = 135 MeV 
and (f3^ pp ) = 0.318, similar to other analyses which have 
neglected decay contributions [^JlTj]. The larger average 
velocity which results in this case is due to the fact that 
all of the observed curvature must now be accounted for 
by transverse flow. 

The high statistical accuracy and large transverse mass 
coverage of the present ir° measurement reveals the con- 
cave curvature of the tt° spectrum over a large mj range, 
which constrains the parameters significantly. This is 
further demonstrated by studying the local slope at each 
ttit- The local (inverse) slope is given by 

-i f f^Y 1 _±_ ( <*M 

Tlocal ~ ~ (fd^ ) \ ) ' (1) 

The local slope results are plotted in Fig. S. Each 
individual value of Ti oca i has been extracted from 3 ad- 
jacent data points of Fig. |l|. The data are compared to 
the hydrodynamical model best fit results of Fig. |l|, as 
well as fits in which the transverse flow velocities have 



1 All limits given use the data for mr — mo > 2GeV/c 2 
as upper limits only to allow for additional hard-scattering 
contributions. 
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FIG. 3. Unnormalized multiplicity distributions as a func- 
tion of the transverse source velocity for the parameter sets 
given in table |. 

been fixed to larger values comparable to those obtained 
by Refs. [jlOj and NA49 § (sets 2 and 3). The corre- 
sponding fit parameters are given in Table [|. The com- 
parison demonstrates that while the large transverse flow 
velocity fits can provide a reasonable description of the 
data up to transverse masses of about 1 GeV, they sig- 
nificantly overpredict the local slopes at large transverse 
mass. While application of the hydrodynamical model at 
large transverse mass is questionable, the model cannot 
overpredict the measured yield. The observed overpre- 
diction therefore rules out the assumption of large trans- 
verse flow velocities, or points to a deficiency in the model 
assumptions used in these fits. 

The curvature in the 7r° spectrum at large transverse 
mass is a result of the distribution of transverse veloci- 
ties. Although the spectrum is not directly sensitive to 
the spatial distribution of particle emission, within this 
model it is dependent indirectly on the spatial distribu- 
tion due to the assumption that the transverse rapidity 
increases linearly with radius. The large curvature at 
large transverse mass is due to high velocity contribu- 
tions which result from the tail of the assumed gaussian 
density profile |Q . Figure || shows the transverse source 
velocity distributions dN/ d(ix for the different parameter 
sets. ^ The curves labelled 1-3 correspond to the calcu- 
lations in figure || using a gaussian spatial profile. In 
addition, velocity profiles are shown for a uniform den- 
sity profile (set 4) and for a Woods-Saxon distribution: 

P(r) = T- ^ TT-TT (2) 

1 + exp [{r - r )/AJ 



2 More precisely these are source emission functions inte- 
grated over all variables except the transverse velocity and 
the rapidity, i.e. they are weighted with the produced parti- 
cle multiplicity. 



with A/ro = 0.02 (set 5). These are included in figures 
H and ||. It is seen that the uniform density assumption 
truncates the high velocity tail resulting in less curvature 
in the pion spectrum, while the Woods-Saxon has a more 
diffuse edge at high (3t- 

While the gaussian and uniform density assumptions 
have very different velocity profiles, it is interesting that 
both can provide acceptable fits to the pion spectrum 
with best fit results with similar (fir) and T parameters, 
which give similar effective temperatures, and which have 
similar velocity widths, Prms, as shown in Table |. Com- 
pared to the gaussian profile result, the best fit result us- 
ing the uniform profile gives a lower temperature of 178 
MeV and would lead to weaker limits of (P^ pp ) = 0.42 
and T low = 134 MeV. Limits cannot be set using the 
Woods-Saxon profile due to increased fit ambiguity. If 
the data for tut — m > 2GeV/c 2 is used only as up- 
per limits, as explained above, a best fit result with 
T = 129 MeV and (f3 T ) = 0.42 is obtained. 

The data presented here can be well described with 
high thermal freeze-out temperatures, similar to temper- 
atures which have been extracted for chemical freeze-out 
|0|, and small transverse flow velocities. f\ On the other 
hand, if the larger velocities obtained in other analyses 
which have considered limited particle spectra together 
with HBT results (||l(],[l3|] persist, then the present anal- 
ysis suggests much lower thermal freeze-out tempera- 
tures. For example, none of the different velocity profile 
assumptions used in this analysis allowed to reproduce 
the results of ref. || - all profiles studied require a tem- 
perature of 90 MeV or less, if ((3t) — 0.55 is assumed. 
The present data obviously provide important informa- 
tion on the shape of the freeze-out velocity distribution. 
A more extensive systematic study would require further 
guidance from full hydrodynamical calculations, which is 
beyond the scope of this paper. Recent hydrodynamical 
model calculations fl^,^(J haw found reasonable agree- 
ment with transverse mass spectra within a broad range 
of assumptions. However, in these studies it was not at- 
tempted to limit the model parameters or assumptions 
by a rigorous comparison to the data. 

In summary, we have argued that hydrodynamical 
models which attempt to extract the thermal freeze-out 
parameters of relativistic heavy ion collisions must pro- 
vide an accurate description of the pion spectra, since 
pions most directly reflect the thermal evironment in 
the late stage of the collision. In particular, models, 
or parameter sets, which overpredict the observed pion 
yields, even at large transverse mass, can immediately be 
ruled out. We have demonstrated that the high accuracy 
neutral pion spectra with large transverse mass cover- 



3 Note again that chemical and thermal freeze-out are not 
necessarily expected to be the same. 
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TABLE I. Parameters for different hydrodynamical model 
fits to the neutral pion spectrum shown in figures ^ and [| 
The temperature T, average and RMS transverse flow ve- 
locity (Pr) and Prms are given together with the effective 
temperature T eff = 77^/(1 - (#r))/(l + (#r)). 



Set 


spatial 
profile 


T 
(MeV) 


(Pt) 


/3rms 


Teff 

(MeV) 


X 2 /dof 


1 


Gauss 


185 ±4 


0.213 ±0.020 


0.107 


230 


25.9/18 


2 


Gauss 


75 ±1 


0.469 


0.199 


125 


386/19 


3 


Gauss 


49 ±1 


0.527 


0.213 


88 


578/19 


4 


Uniform 


178 ±13 


0.274 ±0.046 


0.093 


235 


33.3/18 


5 


WS 


146 t£ 


U.DOO -0.069 


0.137 


214 


26.7/18 



age can constrain the thermal freeze-out parameters and 
model assumptions. Within the context of the hydro- 
dynamical model of Ref. ||, the default velocity profile 
favors large thermal freeze-out temperatures similar to 
the chemical freeze-out temperature determined for the 
same system [Q. Only special choices of the velocity 
profile allow large average freeze-out velocities similar to 
those extracted from other recent analyses which consider 
also HBT results [p| JIo| , p^ |. On the other hand, the cor- 
responding freeze-out temperatures are then « 90 MeV, 
significantly lower than other estimates. The present re- 
sults indicate that the determination of the freeze-out 
parameters remains an open question, ft will be impor- 
tant to determine whether full hydrodynamical models 
can reproduce the high precision pion data and thereby 
constrain the assumed freeze-out hypersurface. 
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